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Abstract. We formulate Grover’s unstructured search algorithm as a chiral quantum 
walk, where transitioning in one direction has a phase conjugate to transitioning in 
the opposite direction. For small phases, this breaking of time-reversal symmetry is 
too small to significantly affect the evolution: the system still approximately evolves 
in its ground and first excited states, rotating to the marked vertex in time ■k'/N /2. 
Increasing the phase does not change the runtime, but rather changes the support 
for the 2D subspace, so the system evolves in its first and second excited states, 
or its second and third excited states, and so forth. Apart from the critical phases 
corresponding to these transitions in the support, which become more frequent as the 
phase grows, this reveals that our model of quantum search is robust against time- 
reversal symmetry breaking. 


PACS numbers: 03.67.Ac 


1. Introduction 


Continuous-time quantum walks BE] are analogues of continuous-time classical Markov 
chains. In them, probability amplitude flows from one vertex of a graph to another by 
evolving by Schrodinger’s equation 


with Hamiltonian 


^walk 2^-^: 


( 1 ) 


where L = A — D is the graph Laplacian, which is composed of the adjacency matrix 
A for which Aij = 1 when vertices i and j are connected (and zero otherwise) and the 
diagonal degree matrix D for which Da = deg(i). The real parameter 7 is the jumping 
rate, or amplitude per time. This evolution causes probability amplitude to transition 
between adjacent vertices. 

For example, consider the quantum walk on the complete graph of N vertices, an 


example of which is shown in hgure la Then the vertices label computational basis 


states {|1),..., \N)} of an iV-dimensional Hilbert space. Since each vertex is connected 
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Figure 1. (a) The complete graph of iV = 5 vertices, (b) Probability overlaps of |s) 
and |w) with the ground and hrst excited states |'0o,i) of H for quantum walk search 
with N = 1023. (c) Success probability for quantum walk search with N = 1023 and 

l = lc = l/N. 


to the N — 1 others, the Laplacian is a.n N x N matrix with —{N — 1) on the diagonal 
and 1 everywhere else: 

/-(iV-1) 1 ■■■ 1 \ 

1 -(N-1) ■■■ 1 

L = 

V 1 1 -(iv-1) 

If the system starts in the equal superposition over all the vertices 

1 ^ 

1=1 

then the system stays in this equilibrium state, up to a global phase, when evolving by 
the quantum walk. 

To construct a search algorithm from the quantum walk, we include an additional 
term in the Hamiltonian 

H =—'yL — \w){w\, (2) 

which marks a basis state |tc) and acts as an oracle [3]. The goal is to hnd |tc) by evolving 
with this Hamiltonian, and we begin in |s) to express our initial lack of knowledge of 
which vertex |tp) might be. Note that the complete graph is homogeneous, so which 
particular vertex is marked does not change the dynamics. Also, each vertex is connected 
to every other, and there is no additional structure to the problem. Thus this is the 
quantum walk formulation of the unstructured, black-box search problem that Grover’s 
algorithm P] famously solves. 

This formulation of Grover’s algorithm is quite easy to solve, as shown by Ghilds 
and Goldstone |5], and which we summarize here. Since the unmarked vertices evolve 


identically, the system evolves in a 2 D subspace. As shown in hgure IT, when 7 takes 
its critical value of 7 c = l/N, this 2D subspace is spanned by eigenstates of H 
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(a) (b) 

Figure 2. (a) Graphic notation for chiral quantum walks. A dotted arrow from i to j 
denotes a directed edge from i to j with weight e‘®, and from j to i with weight e“‘®. 
(b) Chiral quantum walk on the complete graph with N = 5. 



which have an energy gap of AE = Ei — Eq = 2/\/N. Thus the system evolves from |s) 
to |ta) with probability 1 in time tt/AE = iry/N /2, and a detailed calculation showing 
this is available in Section III of [B]. The evolution of the success probability with time 
is shown in figure Ic, and it reaches 1 at 7rA/l023/2 50.241, as expected. 


Such quantum walks are the basis of several of quantum algorithms [Hi and 
even universal computation [H], and they preserve time-reversal symmetry, meaning the 
probability at each vertex at time t is the same as at time —t. Recently, Zimboras et al. 
introduced chiral quantum walks uni, which break time-reversal symmetry by walking 
on directed, weighted graphs, where transitioning in one direction has some phase e^®, 
and transitioning in the opposite direction has the conjugate phase e“^^. We introduce a 
graphic notation for this in figure where a dotted arrow from vertex i to j indicates 
a directed edge from i to j with weight e'®, and from j to i with weight e“‘®. 

By breaking time-reversal symmetry, Zimboras et al. attained significant speedups 
in quantum transport, indicating that such breaking may be a useful tool in engineering 
hardware to efficiently transport energy or information. While faster transport with 
a directed quantum walk mi is unsurprising since allowing transitions in just one 
direction provides immense control over the flow of amplitude, speedups with chiral 
quantum walks are novel because the conjugate relationship that opposite directions 
must maintain prevents one from eliminating one direction of the walk entirely. Chiral 
quantum walks have also been investigated in discrete-time [12], but here we focus on 
the continuous-time formulation. 

Given the speedups obtained by Zimboras et al. in quantum transport, one might 
hope to achieve similar speedups algorithmically. In this paper, we conduct the first 
algorithmic investigation of chiral quantum walks, showing how breaking time-reversal 
symmetry affects quantum search. While the 0{'/N) optimality [TB] of the continuous¬ 
time search Hamiltonian (|^ prevents time-reversal symmetry breaking from improving 
the runtime scaling, it is possible that it improves the constant factor. Of course, it is 
also possible that it makes the search slower. We show that neither of these occur, that 
the runtime is unchanged from tt^/N /2 for large swathes of 6. In spite of this, the value 
of 7 must be chosen differently, and the system may not evolve in its lowest two energy 
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eigenstates. 


In particular, we modify the complete graph as shown in hgure 2b Assuming 


throughout the paper that N is odd, we arrange the vertices in a circular pattern. Then 
using our graphic notation in hgure for chiral quantum walks, each vertex has a 
dotted arrow going from it to the {N — l)/2 vertices ahead (in a clockwise direction) 
of it, and the {N — l)/2 vertices behind it have dotted arrows coming into the vertex. 
Thus each vertex has half its dotted arrows going out and half coming in, but in such 
a way that the graph is homogeneous, so it makes no difference as to which vertex is 
marked for search. 

To further illuminate the chiral quantum walk we have dehned on the modihed 


/ 


L = 


e [2b| the graph Laplacian L for it when N 

= 5 is 

Ni cos 9 

e 

e-ie gi6» 

ei 0 

eie 

—Ni cos 9 

e-ie Q-iO 

eie 

ei 0 

eie 

—Ni cos 9 e~'^ 

e-ie 

e-i 6 » 

ei 6 » 

e’® —Ni cos 9 

e -10 

e-i 6 » 

e-i 6 » 

^lO gi6» 

—A^i COS0 


\ 


V 

where Ni = N — 1. Generalizing this to arbitrary N is straightforward—the hrst column 
starts with —{N — 1) cos6, is followed by (A^ — l)/2 terms that are e'^, and ends with 
[N — l)/2 terms that are e“‘^. The remaining columns are given by cyclic permutations 
of this hrst column, so L is a circulant matrix [TT]. Note that the equal superposition 
|s) is still the equilibrium distribution of this chiral quantum walk ([^, as in the case of 
the normal {6 = 0 ) quantum walk. 

The symmetry of this graph Laplacian L restricts the values of 6 we must consider. 
For example, when 9 > n, then e^^® = so we can simply use 9 — tt < tv and 

reverse all the arrows in hgure which is an equivalent dehnition of the walk. Thus we 
can restrict 9 to the interval [0,7r]. Furthermore, when 9 > 7r/2, then e^‘® = — 
so we can instead use tt — 0 < 7r/2, reverse the arrows in hgure and hip the sign of 
7 in ([^. Thus we can restrict 9 to the interval [0,7r/2], which we do for the rest of the 
paper. 

To turn this chiral quantum walk into a search algorithm, we again include an 
additional oracle term in the Hamiltonian ([^ that marks the vertex to search for. When 
9 > 0, the unmarked vertices no longer evolve identically, and so the system no longer 
evolves in an exact 2D subspace. Thus hnding the runtime and success probability of the 
algorithm by hnding the eigenstates and eigenvalues of H is non-trivial. Nonetheless, 
we solve this problem using the method of Childs and Goldstone [5] , where various sums 
involving the eigenvalues of the negative Laplacian —L can be used to determine the 
critical 7 , as well as the runtime and success probability of the algorithm. We summarize 
this method in the next section. Afterward, we solve search by chiral quantum walk, 
showing that the system still evolves in an approximate 2D subspace, but the value of 9 
can cause the support of the subspace by the eigenvectors of H to change. Away from 
the critical 0 ’s at which the support changes, breaking time-reversal symmetry does not 
affect the search algorithm, so it still succeeds with probability 1 at time ny/N / 2 . 
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We begin by reviewing Childs and Goldstone’s [S] method, which they used to solve 
spatial search on arbitrary-dimensional cubic lattices by quantum walk. We will not 
rederive it, but simply state how to use it. Recall that the search Hamiltonian (|^ has 
two terms: the quantum walk term — 7 L ([^ and the oracle — |tc)(tc|. In general, the 
eigenvalues and eigenvectors of the search Hamiltonian (|^ are hard to hud, but often 
they are easy to hnd for the walk term alone because of the symmetry of the graph. 

In particular, let £j be the eigenvalues of —L, with j = 0,1,..., iV — 1 labeling the 
eigenvalues. Then consider the sum from Childs and Goldstone’s (33): 


5. 


1 

N 


E 

£j^0 


1 

w 


( 3 ) 


From their analysis in Section IV.B., the critical 7 corresponds to Si for large N: 


7c ~ ^ 1 . 


( 4 ) 


As shown in hgure lb, when 7 < 7c, the initial equal superposition state |s) is 
approximately an eigenstate of H, and when 7 > 7c, it is approximately another 
eigenstate of H. Thus the overlap of |s) with the eigenstates of H exhibits a phase 
transition at 7 c. 

As an example, for the normal quantum walk on the complete graph (which is the 

chiral quantum walk with 6 * = 0), —L has eigenvalues 0 with multiplicity 1 and N with 

multiplicity — 1 [15] . Then 

1 N-1 1 

_ 


Si = 


N N 

which is the expected 7 c from the introduction. 

From Childs and Goldstone’s (96), the success probability is approximately 


P* 


at time 


~ - 


vr a/S^ 


2 


Vn. 


( 5 ) 


( 6 ) 


Continuing the above example of normal quantum walk search on the complete graph, 

we already have Si. The other sum is 

1 N-1 1 

_ 


So = 


N m 

Thus 81 / 1/82 ~ 1, and so the success probability roughly reaches 1 at time Try/N /2, as 
expected from the introduction. 


3. Chiral Quantum Search 

Now we apply the method of Childs and Goldstone, which we summarized in the previous 
section, to the chiral quantum walk search algorithm. 
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To calculate the sums Si and S' 2 , we need the eigenvalues of —L. Since —L is a circulant 
matrix, its eigenvalues Sj can be expressed in terms of roots of unity ojj = P|: 

e, = (N-l) cos B - 


where j = 0 , 1 ,..., — 1 , 


Sj — ojjEoJj + ... + oj j ^ 

is the sum of the roots of unity in the upper-half complex plane, and 

N-\ I 1 iV-1 , o 

is its conjugate and corresponds to the sum of the roots of unity in the lower-half complex 
plane. Note that loq = 1, and so Sq = 0- Since Ej is a geometric series, when j ^ 0, 




bJl 


UJc 


-+1 


OJi 


To further simplify this sum, we multiply the top and bottom of it by 1 — cj*: 


Ui 


S,= 


a;/ -c/ 


-+i 


- 1 


Uj -|- 2 i Im (^ 1 — 1 


2 — {uj + u*) 2 — 2 Re{ujj] 

1 sm(gy)+2sm(gyjyi) l a, 

2 2 - 2 cos(«) 2 2 ’ 


where we have dehned 

sm(gy)+2sm(gyjyi) 

' l-cos(M) 

We will further simplify aj in a moment, but hrst let us show how this expression for 
Ej simplihes our expression for the eigenvalues of —L with j 7 ^ 0: 

£,^„ = (/V - 1) cose - e-“ (-1 + - e'" “ if) 

= (TV - 1) COS 0 + ^ (e'® + e-‘^) + i^ (e'^ - e"''') 

= (A^ — 1 ) cos 9 + cos 9 — aj sin 9 

= N cos9 — ajsm.9. (7) 


Thus we have a simple expression for the eigenvalues of — L. 

Now let us return to simplifying aj. For its hrst sine term, we use sin(2M) = 
2sin(M) cos(n). For its second sine term, note that it equals sin( 7 rj — nj/N), so we can 
use sin(M -f- n) = sin(M) cos(n) — sin(n) cos(m), for which the hrst piece is zero: 

_ 2 sin (y) cos (f) - 2 sin ( = ) cos(irj) 

1 - cos (^) 
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Figure 3. The sum (solid black) and its approximation l/(A^cosd) (dashed red) 
as a function of 9 with N = 1023. 


Now we factor 2 sin( 7 rj/A^) in the numerator and use 2sin^(M/2) = 1 — cos(m) in the 
denominator: 


^ 2 sin (f) [cos (f) - cos( 7 rj)] ^ cos (f) - cos( 7 rj) 

2 sin 2 (f) sin(f) 

Note that cos( 7 rj) = 1 when j is even and —1 when j is odd. So we write each case 
separately: 


cos 


*^j,odd 


a 


(¥)-l 


sm 

cos 


j,even 


sm 


ui] 

N j 


N 


= cotl^ 


N 


N 


CSC 


Plugging these into (j^ gives simple expressions for £ 


WO- 


( 8 ) 

( 9 ) 


3.2. Critical Thetas 


Now that we have the eigenvalues £j of —L, we can hnd the sum ([^ 




which is approximately the critical 7 for large N (|^. Since £q = 0, we exclude the j = 0 
term in the sum. But for certain “critical” values of 9, there are other eigenvalues £j^Q 
that are also zero, at which Si diverges. This is shown in hgure Si diverges at several 
critical values of 0, and they are more frequent as 9 approaches 7 r/ 2 . 

We can analytically hnd the 9£s corresponding to these divergences. Since aj^odd 
(|^ is always positive and aj^even is always negative, £j^o ([^ can only be zero when 
j is odd. Setting Sj^odd = 0 and solving for 9, we get 


9c = tan ^ [- J . 

V ®j,odd / 

This can be further simplihed when j scales less than N by Taylor expanding aj^odd ~ 
2N/{'Kj) for large N. Then the critical ^’s are 

9c ^ tan-i 1.0039,1.3617,1.4442,1.4801,1.5002,.... 
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Since j is odd in this expression, let us call these values 6 'ci, 6 'c 3 , 6 'c 5 , and so on. Note that 
they are independent of iV, and they correspond to the divergences of Si in hgure 


3.3. Critical Gamma 


Now let us return to hnding the sum Si, which is approximately the critical 7 for large N 
Q. We assume that 6 is away from its critical values. Otherwise, Si would diverge, and 
7 c would also diverge, which is unphysical since 7 is the jumping rate of the randomly 
walking quantum particle. In particular, we prove below for large N that away from 
the ^c’s. 




( 10 ) 


N cos 9 ’ 

as shown in hgure Note when 0 = 0, this yields 7 c = 1/iV, as expected from the 
introduction. 


using Q, we can rewrite the sum as 
1 

N 


To begin the proof of (10), hrst note from ([^ and (|9j) that a^-j = —C(j- Then 

1 


7=1 ^ 

J 2 

The summand takes the form 

1 1 2a 


+ 


N cos9 — ajsm.9 cos 0 + a,-sin 0 


a + h a — h a? — 


so the sum is 

^ cos^ 6 — sin^ 0 

■ 1 JV-l J 

J 2 

Let us consider this for different ranges of 0. When it is less than its hrst critical value 
(he., 0 < 0ci ~ 1.0039), then the denominator is dominated by A^^cos^0, yielding 


Si ^2 cos 0 ^ 


2 cos 0 iV — 1 


cos^ 0 cos^ 0 2 


N cos 0 ’ 




which is (10). Now consider when 0 is between two consecutive critical values 0cj and 
0 c(j+ 2 ) with J odd. Note that J is small compared to N, otherwise 0 would approach 
7 r /2 where Si diverges (see hgure |^. In this case, we split Si into two sums with odd 
and even j’s, 

*5*1 odd 2 COS 0 ^ ^ Dt 2 ^ 2 n 

^ cos^ 9 — aj sm 0 
‘S'l^even 2 COS 0 ^ ^ 




cos^ 9 — sin^ 0 ’ 


and show that each contributes l/(2iVcos0) to the total sum. For concreteness, we 
have assumed (N — l)/2 to be odd, but our proof still holds when it is even. Beginning 
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with the odd sum, when j < J, the sin^ 6 term dominates the denominator, while for 
j > J, the N'^ cos^ 9 term dominates: 




odd 


2 cos^ 


E 


-1 


a) sin^ 9 

0=1, 3 ,..., J J 


+ 


E 


cos^ 9 

i=J+2,J+4,...,^ 

We want to show that the second sum in the parenthesis dominates the first. To do 
this, we bound each of them. The first can be upper bounded by replacing 9 by 9ci- 

-1 -1 -vrj 


E 


q; 2 sin^ 9 

j = l,3,...,J J 


< 


E 


a] sin^ 6*, 
7r2 + 4(J+l)2 


cl 


E 


j=l,3,...,J 


51^2 Q 


cl 


IbTT 


Ar2 


where we used aj ~ 2N/{Tij) for small j, sin6*ci ~ t: + 4, and 1 + 3 + ... + J = 
(J + 1)^/4. The second sum in the parenthesis can be lower bounded by replacing 9 by 

^^ci: 

1 1 


E 


j—J +2, J+4,..., - 


cos^ 9 


> 


E 


j=J+2,J+4,. 

TT^ + 4 AT • 


cos^ 9ci 


2J- 1 


16 iV2 

where we used cos6'ci ~ + 4 and that the sum has (iV — 2J — l)/4 terms. From 

these, we see that the second sum in the parenthesis dominates over the first when 
J = o{-\fN), which is true so that 9 is not near its critical values. Thus the leading- 
order behavior of the odd sum is 

1 2 cos 9 N — 2J — 1 


Si 


odd 


2 cos^ 


E 


j—J+2,J+A,..., 


A^2 cos2 9 A^2 cos2 9 


2N cos 9 

Now to complete the proof that Si ~ l/(A^cos6'), we note that the even sum S'!,even is 
always dominated by cos^ 9 in its denominator when 9 does not approach 7r/2. Then 
it is also approximately l/(2iVcos6') for large iV, which when combined with the odd 


sum yields (10). □ 


For nonzero values of 6*, some overlap plots similar to figure [Tb] are shown in figure]^ 


Let us discuss each of them. In figure 4a, we have 9 = 0.8, and our formula (10) yields 
the critical 7 ciV SiN 1.44. This corresponds to where the ground and first excited 
states of the Hamiltonian are proportional to |s) =F |tc), so the system rotates from |s) 


to Iw), as desired. In the next plot, figure 4b shows a similar behavior when 9 = 1.2, 
but now at 7 ciV SiN ^ 2.76 and with the first and second excited states. So we 
still expect the system to rotate from |s) to |w), but the evolution approximately occurs 
in the 2D subspace spanned by I'i/'i) and 1-02) instead of l-i/lo) and li/li)- Similarly in 


figure 4c, for which 9 = 1.4, the 2D support has shifted again, this time to \'tp 2 ) and 
\'ipQ at the critical 7 cA^ ~ S'liV 5.88. Next, we will explain why these shifts occur as 
9 increases and derive the critical 6*’s at which these transitions occur. 





















Quantum Walk Search with Time-Reversal Symmetry Breaking 


10 



Figure 4. Probability overlaps of |s) and |w) with the eigenenergies |'0i) of H for 
search by chiral quantum walk with N = 1023 and (a) 9 = 0.8, (b) 0 = 1.2, and (c) 
0 = 1.4. 



Figure 5. F{E) for = 5, 0 = 0.6, and 7 = 1 . As indicated by circles, F[E) = 1 
corresponds to eigenenergies of the search Hamiltonian. 


3.4- 2D Subspace Support 

To explain why the energy levels that snpport the system change with 6 , we must look at 
some details of Childs and Goldstone’s |5] method, which we ignored in our summary in 
Section]^ From their (26), the eigenvalues Ej of —L can be used to hnd the eigenenergies 
Ea of the search Hamiltonian ([^ by solving 

F{E,) = 1 , 

where 

j •’ 

An example of E{E) is plotted in hgure and we give some properties of it that 
were identihed by Childs and Goldstone. EiET) has poles where E = ^E^^ F'{E) > 0 
everywhere (except at the poles), and E{E) —)■ 0 as —)■ ±cx). So there exists an 
eigenvalue Ea of H between every adjacent pair of poles, and one left of the smallest 
pole. If 9 is sufficiently small so that = 0 is the lowest eigenvalue of —L, then the 
ground state energy is negative, and the other eigenenergies are positive. 
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Figure 6. The lowest six eigenenergies of the chiral quantum walk search Hamiltonian 
([^ with N = 1023 and 7 = 


As 6 increases, however, then Q and (|^ imply that 8j decreases when j is odd. 
In particular, £j = 0 at the critical ^’s from earlier: 

9c = tan“^ ( ) 

V ttj,odd / 

tan-i ^ 1.0039,1.3617,1.4442,1.4801,1.5002,..., 

at which a level crossing (degeneracy) occurs in the eigenenergies of H. As 6 continues 
to increase beyond critical values, eigenenergies of H become negative, and the system 
approximately evolves in higher and higher energy eigenstates of H. 

The lowest six eigenenergies of H are shown in hgure as 0 is increased. When 6 is 
sufficiently below the first critical value of 1.0039, the system remains in its ground and 
hrst excited states. Since the energy gap is constant throughout this region, the runtime 
is unaffected by 6. When 6 crosses the hrst critical value so that it is sufficiently between 
1.0039 and 1.3617, the system transitions to being in the hrst and second excited states, 
with the same energy gap as before. So the runtime is unchanged, even though higher 
energy eigenstates of H support the evolution. This continues as each 6c is crossed, until 
the critical values occur so close together that there is no opportunity for the system 
to settle into a 2D subspace. Note that the sharp behavior near each 9c is because we 
used 7 = Si, which diverges at the critical 6 *’s (recall hgure [^. 

Childs and Goldstone’s method follows the two eigenenergies closest to 0, which is 
evident from the derivation of their method [5]. So as the eigenenergies become more 
and more negative, they “drop out,” and the method continues to track the two correct 
energy eigenstates. 

3.5. Success Probability and Runtime 

From hgure and its description in the previous section, we see that the energy gap 
between the two relevant eigenstates is constant away from the critical 0’s. So the 
success probability should reach 1 at time 7r\QN/2. To prove this analytically, we use 
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Figure 7. The sum S 2 (solid black) and its approximation l/(iV^cos^0) (dashed 
red) as a function of 9 with N = 1023. 


(§ and (§, for which we need the snni ([^ 

Finding the behavior of this for 9 away from its critical values follows the same procedure 


as our proof of Fi (10), yielding 


N2 2 0 - ( 11 ) 

IS ^ cos^ 9 

A plot of S 2 with this approximation is shown in figure Note when 0 = 0, this yields 
S 2 = 1/A^^, as expected from the introduction. 

Now that we have expressions for both Si (10) and S 2 10. the success probability 
(|^ of the chiral quantum walk search algorithm reaches 






1 


N cos 8 


= 1 , 


1 


at time (|^ 


^ 7ly/S^ 
2 Si 


Af2 cos2 8 




So even with time-reversal symmetry broken, the runtime and success probability is 
unchanged, so long as we are away from the critical 0’s at which the eigenvectors of H 
that support the evolution change. 

These success probabilities and runtimes are conhrmed in hgure where we plot 
the evolution of the success probability using the same values of 0 that we used in 
hgure 1^ The success probabilities roughly reach 1 at time 7r\/l023/2 50.241. The 

last hgure with 0 = 1.4 is a little short because our approximations for Si and S 2 are 
less accurate as each critical 0 is passed, so bigger N is needed for the probability to 
near 1. 
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Figure 8. Success probability as a function of time for chiral quantum walk search 
with N = 1023, y = Si, and (a) 9 = 0.8, (b) 1.2, and (c) 1.4. 


4. Conclusion 

Time-reversal symmetry breaking in quantum walks was first exploited for quantum 
transport, yielding significant speedups. Here, however, we considered the hrst 
algorithmic investigation of such chiral quantum walks, dehning an unstructured search 
problem on a modified complete graph that retains the homogeneity of the search 
problem. While breaking time-reversal symmetry does affect the jumping rate 7 that 
must be used, and can also change the energy eigenstates in which the system evolves, it 
does not change the runtime of the search for large swathes of 6 . Thus quantum search 
is robust against time-reversal symmetry breaking in the model we have introduced. 

Acknowledgments 

Thanks to Jacob Biamonte for useful discussions. This work was supported by the 
European Union Seventh Framework Programme (FP7/2007-2013) under the QALGO 
(Grant Agreement No. 600700) project, and the ERG Advanced Grant MQC. 

References 

[1] Farhi E and Gutmann S 1998 Phys. Rev. A 58(2) 915-928 

[2] Kempe J 2003 Contemp. Phys. 44 307-327 

[3] Mochon C 2007 Phys. Rev. A 75(4) 042313 

[4] Grover L K 1996 A fast quantum mechanical algorithm for database search Proc. 28th Annual 

ACM Symp. on Theory of Computing STOC ’96 (New York: ACM) pp 212-219 

[5] Childs A M and Goldstone J 2004 Phys. Rev. A 70(2) 022314 

[6] Wong T G 2015 arXtv:1502.04567 [quant-ph] 

[7] Shenvi N, Kempe J and Whaley K B 2003 Phys. Rev. A 67(5) 052307 

[8] Farhi E, Goldstone J and Gutmann S 2008 Theory Comput. 4 169-190 

[9] Childs A M 2009 Phys. Rev. Lett. 102(18) 180501 

[10] Zimboras Z, Faccin M, Kadar Z, Whitfield J D, Lanyon B P and Biamonte J 2013 Sci. Rep. 3 

2361 

[11] Hoyer S and Meyer D A 2009 Phys. Rev. A 79(2) 024307 






Quantum Walk Search with Time-Reversal Symmetry Breaking 


14 


[12] Lu D, Biamonte J D, Li J, Li H, Johnson T H, Bergholm V, Faccin M, Zimboras Z, Laflamme R, 

Baugh J and Lloyd S 2014 arXiv:1405.6209 

[13] Farhi E and Gutmann S 1998 Phys. Rev. A 57 2403-2406 

[14] Davis P 1994 Circulant Matrices 2nd ed Chelsea Publishing Series (AMS Chelsea Publishing) 

ISBN 9780828403382 

[15] Chung F R K 1997 Spectral Graph Theory [CBMS Regional Conference Series in Mathematics 

no 92) (American Mathematical Society) 



